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Abstract. We study a scenario of the ellipsoidal universe in the brane world 
cosmology with a cosmological constant in the bulk . From the five-dimensional 
Einstein equations we derive the evolution equations for the eccentricity and the scale 
factor of the universe, which are coupled to each other. It is found that if the anisotropy 
of our universe is originated from a uniform magnetic field inside the brane, the 
eccentricity decays faster in the bulk in comparison with a four-dimensional ellipsoidal 
universe. We also investigate the ellipsoidal universe in the brane-induced gravity and 
find the evolution equation for the eccentricity which has a contribution determined 
by the four- and five-dimensional Newton's constants. The role of the eccentricity is 
discussed in explaining the quadrupole problem of the cosmic microwave background. 



PACS numbers: 04.50.+h, 98.80.Cq 
1. Introduction 



The WMAP three-year results show that the CMB anisotropy data are in a remarkable 
agreement with the simplest inflation model, but interestingly the large-scale feature 
still warrant further attention p]. The suppression of power spectrum at large angular 
scales (8 > 60°), which is reflected in the most distinguishable way in the reduction of the 
quadrupole C2, remains unexplained by the standard inflation model. Several authors 
suggest that the low multipole anomalies in the CMB fluctuations maybe a signal of 
a nontrivial cosmic topology [21 [3j H]. More precise measurements of WMAP showed 
that the quadrupole C2 and octupole C3 are unusually aligned and are concentrated in 
a plane inclined about 30° to the Galactic plane [5J. This motivated an asymmetric 
expansion universe model, in which one direction expands differently from the other 
two (transverse) directions of the equatorial plane [6]. It was further found that if the 
large-scale spatial geometry of our universe is plane symmetric with an eccentricity at 
decoupling of order 10~ 2 , the quadrupole amplitude can be drastically reduced without 
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affecting higher multipoles of the angular power spectrum of the temperature anisotropy 

m. 

In this paper, we explore a scenario of the ellipsoidal universe in the brane world. 
The brane cosmology of a 3-brane universe in a five dimensional spacetime has been 
investigated in Refs. [HI M, HDJ- The Friedmann equation for such a brane cosmology 
shows that the square of the Hubble parameter H depends quadratically on the brane 
energy density, whenever it depends linearly on the matter energy density in the 
standard cosmology [IP] . We now extend the isotropic 4D spacetime to an anisotropic 
spacetime in the brane cosmology, find the equations for the isotropic scale factor as well 
as the eccentricity for anisotropy and then relate the eccentricity with the quadrupole 
anisotropy of CMB. 

The organization of this paper is as follows. In Sec. II, we derive the basic 
equations of a 4D ellipsoidal universe and find that for the isotropic part of the 
total energy-momentum tensor, the evolution of energy density is described by p 1 oc 
(1 — e 2 ) _ ^ +1 ^a~ 3 ^ +1 \ where e is the eccentricity of the universe and a; is a parameter 
for the equation of state. In Sec. Ill, we find the 5D Einstein equations for an ellipsoidal 
universe with a cosmological constant in the bulk. And we obtain the equations 
governing the evolution of the scale factor and the eccentricity and then discuss their 
cosmological consequences. In particular, we show that if anisotropy of the our universe 
is mainly originated from a uniform magnetic field inside the brane, the eccentricity 
decays faster in the bulk. In Sec. IV, we present the evolution of the eccentricity in the 
brane-induced gravity and find some higher order contributions due to the presence of 
an extra dimension. In Sec. V, we discuss the relation between the eccentricity in the 
brane and the CMB anisotropy and conclude the paper. 

2. Four dimensional planar symmetric universe 

In this section, we derive the basic equations that describe the evolution of a 4D 
ellipsoidal universe. The general plane-symmetric metric is given by [11] 

ds 2 = dt 2 - a 2 (t)(dx\ + dxj) - c 2 (t) dx\, (1) 

where the scale factors a and c are functions of the cosmic time t only and x\x^ denote 
the coordinates of the plane of symmetry. In terms of the eccentricity defined by 

the metric can be rewritten as 

ds 2 = dt 2 - a 2 (t)(dxl + dx\) - (1 - e 2 (t))a 2 (t)dx 2 3 . (3) 

Notice that we have tentatively assumed that a > c. In fact, whether the shape of 
the universe is an oblate (a > c) or prolate (a < c) spheroid depends on the form of 
anisotropic energy-momentum density. The eccentricity for a prolate sphere (a < c) is 

given by another form, e = y 1 — (f) 2 - Here, we should justify the use of an oblate 
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Figure 1. A uniform magnetic field or a cosmic string results in an oblate spheroid 
shape of the universe (left figure), while domain walls lead to prolate spheroids (right 
figure) [6]. 



spheroid and the corresponding eccentricity of the form ([2]). As we have mentioned the 
universe would have expanded isotropically before an anisotropic expansion would have 
become important. The isotropic tension densities plus anisotropic tension densities 
would have caused the spherically symmetric sphere evolve into a spheroid. At a later 
stage of the evolution of the universe, when all of the contributions except the vacuum 
energy (cosmological constant) faded away, the longitudinal and transverse directions 
expand in an equal ratio and the expansion became isotropic. It was proved in Ref. [6] 
that if the isotropy symmetry of the universe is broken by a uniform magnetic field or a 
cosmic string, then the resulting shape of the universe is an oblate sphere (see figure 1). 
In our discussions below, we focus mainly on a homogenous but anisotropic universe, 
where the anisotropy is contributed by a uniform magnetic field. 

The energy-momentum tensor of the whole universe can be in general given by 

T% = diag (p, -p\\,-p\\,-p±). (4) 

Then the Einstein equations read 

a\ 2 2a ee 8n ^ ,_ N 
a J Sal — e l 6 

a\ 2 2a h ee ee e 2 _ , . 

a) + ^- 3 aT^-—*-a^ = - 8nG ^ (6) 

-) 2 + 2- = -^G P±1 (7) 
a/ a 

where overdots denote derivatives with respect to the cosmic time. ^From Eqs. (J5J), fl§D 
and (j7j), we get the conservation equation for the energy- momentum tensor 

e + 2Q( P+P ,) + (i-2 T ^)( P + P± )=o. (8) 
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In principle, the total energy-momentum tensor, T^, can be separated into two different 
parts: an anisotropic contribution (T^)^ = diag(p A , — p A , — p A , — p A ), which may 
include magnetic fields or static aligned strings or static stacked walls, and an isotropic 
contribution, (T/)^ = diag^p 1 , —p 1 , —p 1 , — p 1 ), which includes symmetric contributions 
from vacuum energy or matter or radiation. As in the isotropic universe, we may have 
a thermodynamic relation 

(9) 

where T is the temperature. Up to an additive constant, the entropy in a volume V is 
given by 

s = (p i + p 1 )^. (10) 



Using V oc y/1 — e 2 a 3 , we obtain 

| = (1 _ e 2 )-t(-2ee)a 3 + 3a 2 a(l - e 2 )^ + a 3 (l - e 2 )\ (-f^f)- (U) 

Assuming an adiabatic expansion of the universe, where the entropy in a comoving 
volume is conserved, from Eq. ffTTl) we obtain a conservation equation 

P 1 + 3 (-) (/ + /) - ^~ 2 {P 1 + P 1 ) = 0. (12) 
\a/ 1 — e z 

In the limiting case of e = 0, Eq. (TT2|) is exactly the well-known conservation equation 

for the energy-momentum tensor in the standard cosmology. For a matter with the 

equation of state p 1 = up 1 , by solving Eq. (fT2l) we find that the density evolves as 

p'cx(l-e 2 )-^ +1 V 3 ^ +1 ). (13) 

Therefore, we have p 1 oc (1 — e 2 ) _1 a -3 in the matter-dominant era (p 1 = 0) and 
p 1 oc (1 — e 2 ) -4 / 3 a~ 4 in the radiation-dominant era (p 1 = p 7 /3). Subtracting Eq. (fT2l for 
the isotropic part from the conservation equation (jSj) for the total energy-momentum, 
we obtain a similar equation for the anisotropic part 

p A + 2 (£) (p^ 4 + pjf ) + - J^L) + ) = 0. (14) 

Some exact solutions of the Einstein equations for kinds of plane symmetric plus isotropic 
components are given in Ref. [6]. 

We now give an example of a uniform magnetic field here, the simplest case, whose 
physical interpretation in the brane-world scenario will be discussed in Sec. IV. Let us 
consider the ellipsoidal universe in the matter-dominant era (p 1 = 0), and normalize the 
scale factors such that a (to) = c(t ) = 1 and thus e(to) = in the present time. The 
uniform magnetic field has the energy-momentum tensor (Ta)^ v = p A diag(l, — 1, —1, 1), 
where p A = B 2 /8tt is the magnetic energy density. Here, it is assumed that the magnetic 
field is frozen into the plasma due to the high conductivity of the primordial plasma 
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and the magnetic field evolves as B oc a 2 . Then, for a small eccentricity and thus 
(1 — e 2 ) ~ 1, Eqs. (jSJ), (JTj) and (TH1) can be approximately written as 

— (ee) + 3iiee = 16nGp A , (15) 
at 

if 2 + 2- = 8vrGp A , (16) 
a 

/ + 4ifp A = 0, (17) 

where if = a/a is the Hubble constant and A is the cosmological constant. In the matter- 
dominant era, we have a(t) ~ (3iiot/2) 2 / 3 and then if = 2/31 The solution of Eq. (fi"5l) 
is e 2 = 8fif 0) (l - 3a- 1 + 2a" 3 / 2 ), where fif 0) = p A (t )/p'&\ where p£ } = 3H 2 /8irG is 
the actual critical energy density. Since a(t) < a(to) = 1, the dominant term of the 
eccentricity is e 2 ~ 160^ a -3 ' 2 . The result of Ref. [7] shows that a small eccentricity 
of order e^ec ~ 10~ 2 at the decoupling epoch generated by the uniform cosmic magnetic 
field with a strength Bq ~ 10~ 9 G can explain the quadrupole problem without affecting 
higher multipoles of the angular power spectrum of the temperature anisotropy. 



3. Ellipsoidal universe in a bulk 

In the following, we discuss the ellipsoidal universe in the brane world cosmology. 
The main purpose of this section is to derive the corresponding equations governing 
the evolution of the ellipsoidal universe for the brane metric. In fact, recent results 
show that brane world cosmological models have the capability to endow dark energy 
with an excitingly new possibility {u < —1) without suffering from the problems faced 
by phantom energy [12]. We consider a 5D spacetime metric with an induced plane- 
symmetric metric on the brane 

ds 2 = gABdx A dx B = g^dx^dx" — b 2 dy 2 , (18) 

where y is the coordinate in the fifth dimension. Here and hereafter, we focus our 
attention on the hypersurface defined by y = 0, which we identify with the world 
volume of the brane that forms our universe. The upper case Latin letters A,B,... 
denote 5-dimensional indices, Greek letters p,, u, ... denote indices parallel to the brane 
world volume, 5 an index transverse to the brane, and Latin letters i,j, ... denote space- 
like indices parallel to the brane world volume. We can further take a metric of the 
form 

ds 2 = n 2 (r, y)dr 2 — a 2 (r, y)(dx\ + dx\) — (l — e 2 (r, y)) a 2 (r, y)dx\ 

-b 2 (r,y)dy 2 , (19) 

where the X1X2 plane is the plane of symmetry and n(r, y) is a lapse function. The 
five-dimensional Einstein equations take the usual form 

Gab = Rab — -^Rqab = k 2 Tab, (20) 

where Rab is the five- dimensional Ricci tensor, R = q ab Rab is the scalar curvature, 
and the constant k is related to the five-dimensional Newton's constant G/g\ and the 
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five-dimensional reduced Planck mass M( 5 ) by the relation, h? = 8nG^ = M, 5 y The 
energy-momentum tensor can be decomposed into two parts 

T b = T B \bulk ~\~T B \brane-, (21) 

where T A B \b u ik is the energy- momentum tensor of the bulk matter. The bulk tensor 
T A B \buik can be further decomposed into two parts: the isotropic contribution, 
{Ti) A B \buik= p^diag(— 1, — 1, — 1, — 1, — 1), from a cosmological constant, and the 
anisotropic part, (Tj) A B \buik= diag(p B ,—P B p—P B ^,—P B± ,—P BT ), where the bulk 
energy density and pressures are independent of the coordinate y. The second term 
T A B \brane corresponds to the matter on the brane(y = 0). The most general energy- 
momentum tensor consistent with the planar symmetry takes the form 

T A B \brane= ^y-diag(p, -p\\,-p\\,-p±,0). (22) 



Substituting Eq. (119]) into the Einstein equation, we get the non-vanishing 
components of the Einstein tensor G A B . 

~ 3 / a' 2 a'b' 3ea'e' a" 
1 ~ P {—^2 a (i _ e 2) ~ — 

a fa b\ eed 
anr \a oJ an z {± — e z ) 
~ x 2 If a' 2 2a'b' 2ea'e' „a!n' b'n' 2a" n" 

G I 1 - G I 2 - To ( -—7 + — T + 



ah ' a(l — e 2 ) an bn a n 
a 2 2ab 2ede „ah 



n 2 \a 2 ab a(l — e 2 ) an 

£ 3 1 / a' 2 2a' b' a'ri b'n' 2a" 
73 b 2 \ a 2 ab an nb a 

_l_ ^ ^ 2 _|_ 2 _|_ 2 ^ _|_ ^ 1 (25) 

n 2 \a? ab an bn a b 



2 a'n> 




b'n' 




+ 


bn 


an 


bh 


2d 


b 


bn 


a 


+ b 








n"^ 


) 




n , 







(24) 



^ 1 /3n'd ^ 3a'b ea'e ^ 3ee'd 3d'\ 
5 n 2 \ an ab a(l — e 2 ) a(l — e 2 ) a / 
~ 5 1 / 3n'd 3a'& ea'e 3ee'd 3d'\ , . 

I0= P\~i^ ~ab + a(l - e 2 ) ~ a(l - e 2 ) + ~~a~ ) ' ( ^ 

^ 5 1 /3d 2 3eed 3dn 3a\ 
15 n 2 \ a 2 a(l — e 2 ) an a / 
1 / 3a' 2 ea'e' 3a'n'\ 
+ V V "a 2 " + a(l -e 2 ) ~ ~^rJ ' (28 ^ 
In the above expressions, primes stand for derivatives with respect to y, and overdots 
for derivatives with respect to the cosmic time r. The reason why we can separate 
the Einstein tensors into isotropic parts and anisotropic parts is that the corresponding 
parts of the energy-momentum tensor obey their own conservation equations, Eqs. (fT2l 
and (JHj). We assume that there is no flow of matter along the fifth dimension, i. e. 
(Ti)o5 = 0, which in turn implies that G° T 5 = G\ = 0. Then the components (0, 0) 
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and (5, 5) of Einstein equations (I23p and (1281) in the bulk can be written in the simple 
form 

f'~(l-eVf; 0) ^y(l-« 5) (29) 

where F is a function of r and y defined by 

„. . . r/da\ 2 /a'a\ 2 l , . 

F(r, „) = <!- e') [(-) -(-)]. (30) 

Integrating Eq. (!30l) . we obtain the relation 

dx 2 _/a'x 2 K 2 |aXl-e 2 ) f0 d 

aJ "Lov 1 + 6 /0 6a 4 (l-e 2 ) /0+ a 4 (l-e 2 )' [6 ' 

where C\ is a constant of integration. The above equation shows that the scale factor 
a(r, y) and the eccentricity are coupled to each other. If e = 0, then Eq. ( !3~Tj) reduces 
exactly to the result of Ref . [10] . 

We now consider the anisotropic contribution of the total energy-momentum tensor. 
The corresponding anisotropic part of the Einstein tensors can be written as 

1 / „ e' 2 eb'e' ea'e' 



G ° A °- b 2 (l - e 2 )V 6 " + (1-e 2 ) ' + 
1 / ede ebe \ 



n 2 (l — e 2 ) V a b J 

~ x _ ~ 2 1 /ea'e' eb'e' e' 2 , , „ 

^ A 6 2 (l-e 2 )V a 6 (1-e 2 ) 

1 /eae ebe e 2 



/eae ebe e een .. \ ,„ . 

(-T + -T + T, 5T + ee > ( 33 ) 



n 2 (l — e 2 ) V 6 6 (1 — e 2 ) n 
^3 = 0, ^ (34) 

1 / ee'd ee'b ea'a e'e en'e 



G A 5 = -Tn o\ T- + + 71 o\ + ee' , (35) 

^ n 2 (l-e 2 )V a 6 a (1-e 2 ) n J 

~r 1 /ee'd ee'6 ea'd e'e en'e .A 

^»=65(r^)(— + — "IT "(1375)+ — " ee )- (36) 

^ 5 1 /ea'e' ^ ee'n' ede e 2 ^ een ee\ 

5 1 — e 2 V ab 2 b 2 an 2 n 2 n 3 n 2 / 

Here we also assume that 5 = 0, that is to say no flow of anisotropic matter 

(magnetic fields or cosmic strings or domain walls) along the fifth dimension and so 

G\ 5 = G b A = 0. Eqs. (|32|) and (l3~Tj) can be written as 

K' ~ 2ee'(l - e 2 )aVf° , K ~ 2ee(l - e 2 )aVTj 5 , (38) 

where is defined by 

aee'\ 2 /aee\ 2 



Integrating Eq. (38) we obtain 



(( \ 2 / ee' \ 2 1 „- n f(l-e 2 )da 2 - n C 2 , N 



l-e 2 )n/ V(l-e 2 )67 2 ^ u a 2 (l - e 2 ) * u a 2 (l - e 2 ) 2 
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where C2 is another integration constant. 

We now take the brane into consideration by using the Israel's junction condition 
|13j . When the coordinate system (18) is chosen, the extrinsic curvature tensor of a 
given hypersurface (for example, the y = surface) is defined by K^ v = d y g^ u /2. The 
Israel's junction condition at y — can be written as [HI [15] 

.2 



{K, u - Kg^} y = -—T^, (41) 



K 

~2 

where K is defined by K = K^g^, is the energy-momentum tensor of the brane. 
We have assumed a brane with Z2 symmetry and the index y means accordingly that 
the value of the component is taken over one side of the brane, namely at y = + . 
The junction condition on the planar symmetric metric background are given by 

bl yP+ y(P|| ~Px), (42) 

k 2 (p\\-P±), (43) 



ao&o 3 3 

e [e'] 

Ml - eg) 



[n'l k 2 



n obo T (2P + 2 P||+P±) (44) 

where the subscript for a, b, and n means that they are taken at y — 0, and 
[Q] = Q(0 + ) ~~ Q(0~) denotes the jump of the function Q across y = 0. From Eqs. 
(42)-(44), we can see that when py = pj_ the boundary conditions reduce to the spherical 
cases discussed in Ref. [10J. Assuming the symmetry y <-> — y for simplicity, the junction 
conditions (42)-(44) can be used to compute a' and e' on the two sides of the brane, and 
by continuity when y — > —y, Eq. (31) and Eq. (40) are reexpressed as (after setting 



no 



a\i_K 2 p I B t k 4 P 2 ^ ^ 2 Ja 4 d(l-e 2 )j t C x 



+ ^ + 36 {P "- p±) ~ MTl - e 2 ) PB + 



) 6 36 36 6(2q(1 — < ^1 "q\ ' — ■ 01 

( e e Y K \ ^2 , « 2 a , [0- ~ e 2 )da 2 A C 2 

= T ( ^» - P ^ + + og(l - eg)^ ^ + ag(l-e^ ' (46) 

The above two equations are the main results of our work. From Eq. (45), we can see 
that when pu = p±, the scale factor is exactly what was found in Ref. [10]. If the third 
term in Eq. (45) is small compared with the other terms, it play a role of perturbation 
and the Friedman equation does not deviate much from that of Ref. [10J. Eq. (46) 
describes the evolution of the eccentricity on the brane world, which shows that the 
evolution of the eccentricity depends on both the brane anisotropic pressures and the 
bulk anisotropic energy density. We should notice that the energy conservation on the 
brane, which was obtained in Sec. II, still works here 

p + 2(-)(p + m)+(-- (p + P± )=0. (47) 



a j v ' 1 ' \ a I — e 2, 

Since the scale factor and the eccentricity are coupled to each other, we do 
not expect to solve Eqs. (45) and (46) exactly. But for a small eccentricity, i.e. 



Ellipsoidal universe in the brane world 



9 



p ~ (1 — e 2 ) ( w+1 )a 3 ( a;+1 ) ps a 3 ( a;+1 ) and when the third term in Eq. (45) can be 
neglected, the time evolution of the scale factor can be given by [10] 

ao(f) = a^K 2 p^{^K 2 p h t 2 + |t) V9 , (48) 

where the energy density on the brane has been decomposed into three parts p = 
p 1 + p A + Pa, q = 3(u +1), a* and p* are constants, and pa is a constant that 
represents an intrinsic tension of the brane. Here the Randall-Sundrum relation, 
K 2 p^/6 + K 2 p A /36 = has been used [16]. Eq. (48) indicates that at a very early 
universe, the cosmology is characterized by a(t) ~ t 1 ^, while at a late time it is 
described by the standard cosmology, a ~ t 2 l q . We are going to solve Eq. (46), by 
considering only the first term and neglecting others. We still assume the anisotropy of 
the our universe is contributed by a uniform magnetic field with the energy-momentum 
tensor (T^)^ = p A diag(l, —1, — 1, 1), and thus from Eq. (47) we find that p A oc a -4 . 
Substituting a ~ t 1 ^ , H = 1/qt, and p A = p A a~ A back to Eq. (46), we obtain 

S + C» (49) 

q — 4 

where C3 is an integration constant. Similarly, at a late time of the universe, a ~ t 2 l q , 
H = 2/qt, and p A = p A a~ A , the eccentricity is expressed as 

9 2K 2 p A t 1 ~l „ . , 

e 2 = ^° fi +C 4 , 50 
g — 8 

where C4 is also another integration constant. In the matter-dominant era (q = 3), we 
find that the late-time evolution of eccentricity in the bulk approximately is e 2 ~ a~ 5 / 2 , 
which decays a little faster than that in 4D universe, while e 2 ~ a -1 for a very early 
universe a ~ t 1 / 9 with q = 3. 



4. Ellipsoidal Universe in the brane-induced gravity model 

In this section, we explore a scenario of the ellipsoidal universe based on the Dvali- 
Gabadadze-Porrati (DGP) model of brane-induced gravity [T7]. In this model, the 
3-brane is embedded in a spacetime with an infinite-size extra dimension. The usual 
gravitational laws is obtained by adding to the action of the brane an Einstein-Hilbert 
term computed with the intrinsic curvature on the brane. Particularly, one recovers a 
standard four-dimensional (4D) Newtonian potential for small distances, whereas gravity 
is in a 5D regime for large distances. The cosmology of this model in the case of a 5D 
bulk was studied by Daffayet, Dvali and Gabadadze. It is shown there that if the 
cosmological model contains a scalar curvature term in the action for the brane, besides 
the brane and bulk cosmological constraints, the presence of the scalar curvature term in 
the brane action can lead to a late-time acceleration of the universe even in the absence 
of any material form of dark energy 0, [9] . 
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Following the framework of Refs. [5J [9J, we consider a 3-brane embedded in a 5D 
spacetime with an intrinsic curvature term induced on the brane and the action of the 
form 



S {5) = J d b x^/-gR + J d b xC m - d A x^R (51) 

The first term in Eq. (51) corresponds to the Einstein-Hilbert action in five dimensions 
for a 5D metric ^^(bulk metric). Similarly, the last term in (51) is the Einstein-Hilbert 
action for the induced metric g c d on the brane, R being its scalar curvature. The induced 
metric g c d is defined as usual from the bulk metric, g c d = d c X A ddX B c/AB, where X A (x c ) 
represents the coordinates of an event on the brane labeled by x c . The second term in 
(51) denotes the matter content. 

The five-dimensional Einstein equations in the brane-induced gravity is given by 

Gab = Rab — -^Rqab = k 2 Sab, (52) 

where the tensor S is the sum of the energy-momentum tensor T of matter and the 
contribution coming from the scalar curvature of the brane. We denote the latter 
contribution U, so 

S A B = f A B + U A B . (53) 

The junction condition in the brane-induced gravity model is replaced by {K^ w — 
Kg^ u }g = — (k 2 /2)S^ u , with S^ u = T^ v — {\/p 2 )U^ V) where p 2 is the four-dimensional 
Newton's constant, p 2 = 87rG( 4 ) = M,~ 2 . The Israel boundary condition is then written 



as 



K 2 K 2 , , K 2 



-P+-r(P\\ + Px) 



a b 3 ' 3 ^ " 3/i 2 ng 



eo[e'\ 



bo(l-e 2 ) 



a 2 e Q de e e e eh e 2 \ 

6 a 2 + a (l - e 2 ) + 1 - e 2 n (l - e 2 ) + (1 - e 2 ) 2 )' l " } 

K 2 {p\\ -p±) 

k 2 ( 3eoae e e eoeh e 2 



p2 n 2 \a (l-e 2 ) + 1-e 2 n (l - e 2 ) + (1 - e 2 ) 2 > (55) 
\n'\ k 2 . . k 2 (a 2 ah 2d\ 2k 2 

e ae e e e eh e 2 \ 

+ 1 3 - — — + ~ 3^ ) • ( 56 ) 

o) 



ao(l-eg) I -el n (l-eg) (1 - e, 



Substituting Eqs. (31), (40), (54) and (55) into the (0,0)-component of the Israel's 
junction condition and assuming p B , p B , C\ and C 2 to be approximately zero, we have 

where H = a/ao, e is the sign of [a'} such that e = ±1 and the lapse function is set to 
be no = 1. One should note that when p 2 / 'k 2 < 2p 2 / 'k 2 <C H, or in terms of the Hubble 
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radius H^ 1 , when H^ 1 <C M? 4 J2ML the 4D ellipsoidal universe is recovered, i.e. one 
can return to Eq. (5). 

We notice that from Eqs. (40) and (55) with vanishing anisotropic bulk energy- 
momentum tensor p"g and vanishing C 2 > it follows that 

(A + (dW) + (i^y ( 3 i + 2 4) = ^ ~ < 58 > 

When 2ep 2 / k 2 in the second term can be neglected, then Eq. (58) exactly equals to 
Eq. (J71) subtracted by Eq. (ED- In the small eccentricity approximation 1 — e 2 ~ 1, we 
solve Eq. (58) in the matter-dominant era with a ~ t 2 ^ 3 and H = 2/31 We normalize 
the scale factor such that a(to) = 1 and e(t ) = in the present time and consider 
the anisotropy on the brane contributed by a uniform magnetic field with the energy- 
momentum tensor (T^)^ = p A diag(l, — 1, — 1, 1). The solution of Eq. (58) can be 
approximately written as 

e 2 (t)~ J\7cp£e- 2e &'t'- 2 dt' 

pto , .,2 zLmS'2 

= jf 4^(1-26^' + ^ + ...)^. (59) 

Eq. (59) indicates that e should be 1, otherwise the eccentricity in the brane does not 
converge when t — > oo. Fortunately, Ref. [Sj shows that the brane cosmology with 
e = 1 can produce a late-time accelerated expansion. When 2/i 2 //? 2 i7 and then 

e -2 ^* ~ e~ 2(i iz u ~ 1 ; thus e 2 (t) ~ t^ 1 ~ a -3 / 2 , which is exactly the dominant term 
of the 4D case. Integrating Eq. (59), the eccentricity on the brane can be written as 

e'W^^ + s!^-^), (60) 

where tt£ 0) = p A (t )/p^\ p { ° ] = 3H 2 /8ttG and r c = K 2 /2p 2 = M ( 2 4) /2M ( 3 5) is the length 
scale for the crossover between 4D gravity and the 5D gravity regimes [TTJ. Differently 
from that of the 4D case, the brane-induced gravity contributes some higher order 
corrections to the eccentricity duo to the presence of an extra dimension. 

5. Discussion and Conclusion 

In the brane cosmology we have found the Friedmann-like equation for the ellipsoidal 
universe in the bulk, which depends on the geometry and matter (both isotropic and 
anisotropic) content of the brane. We also have found that the evolution equation of the 
eccentricity in the bulk depends only on the anisotropic pressures inside the brane and 
the anisotropic energy density in the bulk, except a constant parameter. The evolution 
equation of the eccentricity is coupled to that of the scale factor. As a model calculation, 
we have shown that if only a uniform magnetic field inside the brane contributes to the 
anisotropy of our universe but the anisotropic energy density from the bulk and other 
terms are neglected, the evolution of the eccentricity decays faster than that in a 4D 
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universe. To compare with the 4D ellipsoidal universe, we have considered the ellipsoidal 
universe with a 3-brane embedded in a 5D spacetime with an intrinsic curvature term 
included in the brane action. The results show that the usual ellipsoidal cosmology is 
recovered for Hubble radii smaller than the crossover scale given by r c = M 2 ^/2M^ 
between 4D and 5D gravity. 

We now briefly discuss the relation between eccentricity and cosmic microwave 
background quadrupole problem. For a small eccentricity of the universe in the brane, 
the metric tensor may be written in a perturbation form 

ds 2 = dt 2 - a 2 (t, y)(5 ij - h ij )dx i dx j - b 2 (t, y)dy 2 , (61) 

where represents a metric perturbation, = e 2 ^^. Suppose a photon emitted 
at t = tdec from the last scattering surface travels along a null geodesic and reaches an 
observer at t — to. Let n l be a unit vector along the null ray from the observer to the 
surface. With respect to the observer at the origin, the photon ray is x l (t) = n(t)n l , 
where r/(t) = f dt' /a{t'). The geodesic deviation between an observer at x l and another 
observer at x % + b~x % along the ray at any time t is given by 51 = (—gijSx l Sxi) oc 
a(l — hijU 1 ^ /2) to first order in h^. ^From the proper velocity v = dSl/dt of one 
observer with respect to another and the redshift of the photon's physical frequency 
Sup/ujp = —v, we find the redshift of the comoving frequency 

Su c S(u p a) 1 dh 



uj c u! p a 2 dt 



y-rfn'St. (62) 



As the temperature of radiation with respect to the comoving observer is proportional 
to the frequency (ST oc u c ), we may find the temperature variation from Eq. (|62|) as 

ST 1 f to 1 



, . ... n % n 3 dt = -e dec n 3 , (63) 

where n 3 (8,4>) = cos 9 cos d — sin 9 sin d cos(<ft — y?) and (p are the angles between 
spherical galactic coordinates and the c-axis and a-axis respectively [7]. 

On the other hand, the relative temperature anisotropy ST(9,(f>)/T leads to the 
power spectrum 



STi / 1 1(1 



yEw 2 > ( 64 ) 



(T) V2-K 21 + 

where ai m is the coefficient of spherical harmonics Yi m (9,4>) and (T) = 2.726 ± O.OlOfT 
is the actual average temperature of the CMB radiation. The power spectrum fl64|) fully 
describes all the CMB anisotropy and I = 2 refers to the quadrupole anisotropy. Recent 
WMAP data hints a violation of statistical isotropy on its largest scales and a missing 
power at scales greater than 60° JIJ. Particularly, the observed quadrupole anisotropy 

(ST 2 ) 2 obs ~ 236/iK 2 , (65) 

deviates by order of magnitude from the prediction of standard inflation 

(ST 2 ) 2 ~ 1252/ifT 2 . (66) 
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This anomaly is called the quadrupole problem. Since we have assumed that the 
large-scale spatial geometry of our universe on the brane is plane symmetric with 
a small eccentricity, the observed CMB anisotropy map is a linear superposition 
of two independent contributions 8T = 5Ta + STj, where 5Ta represents the 
temperature fluctuations due to the anisotropic spacetime background, while 5Tj is 
the standard isotropic fluctuation caused by the inflation-led gravitational potential 
at the last scattering surface. Similarly, we may write a im = af m + af m , where 
a tm = (1/2) J ^ Jq e dec n 3^m s ^ n ^ c ^ c ^ ) - Finally we note that while the quadrupole 
anisotropy contributed by the uniform magnetic field is Sa = <5T 2 /(T) = (2/5 v^e^ 
with an eccentricity of order 10~ 2 , the predicted quadrupole anisotropy can be in a 
range of 46.2/iK 2 < (<5T 2 ) 2 < 1001.6 /xK 2 (for detailed calculations, see Ref. [7]). Thus 
the data is in agreement with observations. 

The equation that governs the evolution of eccentricity in a five-dimensional bulk 
is found in Sec. Ill, which is shown to depend on the anisotropic pressures inside 
the brane and the anisotropic energy density in the bulk. If anisotropy of the our 
universe is contributed by a uniform magnetic field inside the brane, and neglecting 
the anisotropic energy density from the bulk and other terms, the evolution of the 
eccentricity decays faster than that in a 4D ellipsoidal universe. Since the evolution of 
eccentricity here is described by a 5-dimensional Newton's constant and is inconvenient 
to be compared with observational data, we have come to focus on the case of ellipsoidal 
universe in the brane-induced gravity in Sec. IV. The perturbational metric takes the 
form of Eq. (61) where we have assumed the fifth dimension to be stabilized under the 
perturbation. This is because the origin of eccentricity is contributed by the magnetic 
field which according to Hofava-Witten picture is confined to a 3-brane [18]. From 
Eq. (60), we can determine the value of e 2 dec . At the decoupling, t = t dec , we have 

— 3/2 

e L - l n M z dic ~ 2lB ^r ~ f ^H' where Zdec ~ 1088 is the redshift at decoupling [T3]. 
For large r c (for example, r c = I^I^'q^o" 1 [5]), we have e 2 ec ~ |^(o)^dec- Fi nan y; we 
get e dec ~ iTgiO -2 /^ 1 !^, where B = B(t ) is defined by p£ = B 2 /8n, and h ~ 0.72. 
Thus, for B ~ (4 - 5) x 10~ 9 G, we have e dec ~ (0.1 - 0.2) x 10~ 2 . Therefore, in the 
brane-induced gravity model, the evolution of eccentricity in the brane can also be used 
to explain the quadrupole problem of the cosmic microwave background. 

In summary, we have considered a version of ellipsoidal universe in the brane 
world scenario. Here, we emphasize that the ellipsoidal universe discussed above 
is undergoing a homogenous but anisotropic expansion. That is to say, before the 
onset of inflation, a typical region of the universe is homogeneous and isotropic but in 
some regions an asymmetric expansion driven by magnetic fields would have stretched 
out the regions and left with an imprint through the end of inflation [6J. The 
ellipsoidal universe has been used to explain the suppression of quadrupole moment, 
but not all the low multipoles since in the ellipsoidal coordinate all the coefficients 
af m with I > 2 vanish. This can be easily found by calculating the formula af m = 
(1/2) J* e^nlY^sinOdOdcj). In the above discussions, we have also investigated the 
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ellipsoidal universe in the DGP model and used a large r c in obtaining the value of 
eccentricity. We notice that up to now the comparison between data and the DGP 
model based upon the expansion history of the universe is still contradictory. Recent 
data of Supernova and CMB suggested that the self-accelerating branch of the DGP 
model is somehow disfavored [201 EU E2, E31 EH [25], [26] , while the most recent analysis 
of the new 'gold' data set of supernovae [27] and the CMB shift parameter suggested 
that a flat universe is completely consistent with the DGP model |28[ 129]. However, as 
the density perturbations of the DGP model may differ from those of ACDM model, it is 
still unclear at present whether the DGP model is marginally or significantly disfavored 
or not [3ni EH [32] . A full examination of this issue is beyond the scope of this paper, 
but from Eq. ([601 we can see that a large enough r c contributes little to the value of 
the eccentricity. 
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